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Abstract 



A quantum mechanical system S is indirectly controlled when the 
control affects an ancillary system A and the evolution of S is modified 
through the interaction with A only. A study of indirect controUabil- 
i^ ■ ity gives a description of the set of states that can be obtained for S 

with this scheme. In this paper, we study the indirect controllability of 
y ■ quantum systems in the finite dimensional case. After discussing the 

jrt I relevant definitions, we give a general necessary condition for control- 

^ i lability in Lie algebraic terms. We present a detailed treatment of the 

O^' case where both systems, S and A, are two-dimensional (qubits). In 

particular, we characterize the dynamical Lie algebra associated with 

S+A, extending previous results, and prove that complete controlla- 

^ . bility of 5'+^ and an appropriate notion of indirect controllability are 

t^^ ' equivalent properties for this system. We also prove several further 

rp^ ■ indirect controllability properties for the system of two qubits, and il- 

^—v ' lustrate the role of the Lie algebraic analysis in the study of reachable 

• . states. 

m 
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^ : 1 INTRODUCTION 



Indirect controllability of a quantum system S refers to the situation where 
two quantum systems S and A interact, and the externally applied control 
is allowed to influence the evolution of A only, while the state of S is of 
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interest. The systems S and A are referred to as the target system and 
accessor system, or ancillary system, respectively. Therefore, the state of 
S is indirectly controlled through the interaction with A. The problem of 
indirect controllability is to analyze to what extent it is possible to modify 
the state of S with this scheme in various situations. 

The question of indirect controllability is of theoretical interest for the 
further development of quantum control and of quantum information, as 
indirect controllability is a signature of the capability of the dynamics to 
generate entanglement [20J. In fact, if there is no entanglement between 
the two systems, S and A, they evolve separately and indirect controlla- 
bility is impossible. This question is also of practical relevance since, in 
many experimental setups, it is much easier to control the accessor system 
A rather than the relevant part S. There might be several reasons for that. 
For instance, we might be able to use only control fields tuned at the res- 
onance frequency of some of the particles which make up a multi-particle 
systems or we might be able to access only systems in certain locations. As 
a physical example of indirect control, protocols of short-distance quantum 
communication have been described using spin chains [HE]- The state to be 
transmitted is put in interaction with one end of the chain, the controlled 
part, and it is reproduced at the other side of the chain, without physical 
transmission. So, control at one end of the chain realizes the task of state 
transfer. It has been proved theoretically [7j that control at the end of the 
chain is sufficient to realize a reliable transfer. This can be seen as indi- 
rect control. The target system S is the last spin at one end of the chain 
whose state we want to manipulate (and make equal to the one at the other 
end) . The rest of the spin system can be seen as the accessor system. Spin 
chain schemes can be implemented through super-conducting devices based 
on Josephson junctions [S^ and are potential tools for the implementation 
of quantum technologies. 

There is a large literature studying the control of quantum systems, 
starting with establishing the connection with geometric control theory 
[HllIlllIBlEI] until more recent studies for closed (see, e.g., [SSI El [MIS]) 
and open (see e.g., [5l[13l[26]) quantum systems. Nevertheless, an algebraic 
characterization of indirect controllability is still lacking. In fact, most in- 
vestigations on indirect controllability have given conditions so that the full 
system S+A is completely controllable, that is, every unitary transformation 
can be performed on it. This implies in particular that every unitary opera- 
tion can be performed on S. In these cases, controllability can be proved by 
verifying the standard Lie algebraic conditions for closed systems, that is, 
calculating the so-called dynamical Lie algebra (see definition in Section [3]), 
which fully characterizes the dynamics of the total system [lOj. However, 
as we shall prove, there are cases where indirect controllability on S, in an 
appropriate sense, can be achieved even without complete controllability of 
the total system, and this motivates a more detailed study of this property. 



The goal of this paper is to provide a general framework for the treatment 
of indirect controllability of quantum systems, with special emphasis on the 
case of two coupled qubits S and A. Because of the introduction of the 
partial trace (see definition in the next section) to describe the dynamics of 
the system S, much of the machinery of Lie transformation groups cannot 
be used in this case. Nevertheless Lie algebraic techniques play a crucial 
role in determining many indirect controllability properties. 

We briefly summarize this work. In Section [2] we define the mathemati- 
cal set up and give the basic notations and mathematical tools. In Section 
|3]we define the main notions for indirect controllability and give a general 
necessary criterion for indirect controllability at the Lie algebra level (The- 
orem [1]) . In Sections H] and [5] we present a detailed treatment of the case of 
two interacting quantum bits, S and A. In particular, in Section[5]we give a 
description of the dynamical Lie algebra associated with this system. This 
extends a result of [IS] which only gives a sufficient condition for complete 
controllability. Based on the characterization of the dynamical Lie algebra, 
in Section [5] we prove various indirect controllability properties for the two- 
qubit system presented in Theorems HI [SJ in Proposition 15.21 and in Example 
15.11 We draw some conclusions and give directions for further research in 
Section El 

2 BASIC FACTS CONCERNING THE DYNAM- 
ICS OF INTERACTING SYSTEMS 

The states of the target and ancillary systems, S and A, are described by 
density operators po and pA, respectively, that is, positive semi-definite, 
trace class operatorq^ with unit trace on Hilbert spaces S and A, associated 
with S and A, respectively [6l|20]. The total system 5 + A is described by 
the density operator pT on the Hilbert space S (d A. The sets of density 
operators on 5, A, and S (^ A are convex, the boundary consisting of pure 
states, that is, projectors of rank one of the form p = |^)(^|, with \ip) in 
S, A, OT S A, respectivelycl The state of a subsystem {ps or pa) can be 
extracted from the state of the total system pT through the partial trace 
operation, ps = TrAipr) and pA = T^^siPT)- We recall that, if we express 
an Hermitian operator p & S <^ A as 



E' 



aj p] <E) pf, (1) 



^A trace class operator is a bounded operator A such that Ti{A^ A)^ is finite. We shall 
deal with the finite dimensional case where this requirement is automatically satisfied. 

^We shall make use of the Dirac notation of quantum mechanics denoting a general 
vector in a Hilbert space (ket) by |-) and the associated dual linear operator (bra) by {■!. 



with Hermitian operators Pj € S and Pj G A, and real coefficients Uj, the 
partial trace operation is a linear map from the space of Hermitian operators 
on 5 (X" ^ to the space of Hermitian operators on S (A) defined as [20] 



TTAip) := Y^ aj p^Tripf), Trs{p) := ^ a, pfTrip^). (2) 

J 3 

The main property of the partial trace which will be used in the following 
is that, for unitary operators F and G acting on S and A respectively, it 
holds 

Tr^ (f ® GpF^ ® gA = FTta{p)fI (3) 



In our discussion we are going to assume finite dimensional systems S 
and A, that is, both S and A are finite dimensional vector spaces. Some 
of the definitions and notions we shall introduce hold for the infinite di- 
mensional case as well or can be naturally extended. We shall denote by 
ns and ha the dimensions of S and A, respectively, and the total system 
S + A has dimension nsriA, so ps, PA and pT are represented by Hermitian, 
positive semi-definite, unit trace matrices of dimensions ns x 71,5, ua x ua, 
nsnA X nsUA, respectively. 

The evolution of the total system is determined by a Hermitian operator 
Ht, called the Hamiltonian, which can be conveniently separated in four 
parts as follows o 

Ht = Hs + Ha + Hi + Hc. (4) 

The Hamiltonian Hs determines the evolution of the system S alone. It is of 
the form Hs = B (^ I^a' ^^^ some Hermitian operator B on S. Analogously 
Ha = Ins '^ ^ ^^^ some Hermitian operator C on A. Hj models the two- 
body interaction between S and A, and it has the general form 

m 

Hj = Y,B,(^C„ (5) 

i=i 

for Hermitian operators Bj and Cj on S and A respectively. Finally the 
Hamiltonian He models the interaction with the external control. Since in 
the indirect controllability scheme we assume only control on A, He will be 
written as 

n 

He = Yuk{t)lns ® Lk, (6) 

fc=l 



^Without loss of generality we shall assume in the following that all Hamiltonians have 
zero trace, since we can always decompose a Hamiltonian H as H — — Tr{H)ln + H where 
H has zero trace and n is the dimension of the underlying system. The term ^ Tr(_ff)l„ 
affects the dynamics only for a global phase factor which is irrelevant in quantum mechan- 
ics, and thus can be neglected. Here and in the following 1 denotes the identity operator. 
We shall write 1„ if the identity operator acts on a space of dimension n, when we want 
to highlight such dimension. 



for some control functions Uk{t) and Hermitian operators Lfc|f| Assuming 
units such that h = 1, the total system S + A evolves according to the 
Schrodinger operator equation 

X = -iHT{u{t))X, X(0) = InsUA, (7) 

and the density operator px varies as 

PTit)=Xit)pT{0)xHt). (8) 

The initial state is assumed to be uncorrelated, that is pr(0) = psi^) ® 
Pa{0)- By combining this evolution with the partial trace ^, we obtain the 
dynamics of the relevant system S, 

ps{t) = Tta [X{t)ps{0) pA(0)Xt(t)) . (9) 

In the following, we shall study the map, for given pA, 

Ps^Tia{Xps^PaX^), (10) 

as X varies in the set of all unitary operators available as solutions of ([7]) 
with suitable control actions. 

3 INDIRECT CONTROLLABILITY FOR GEN- 
ERAL SYSTEMS 

We now give several definitions of controllability, that might be appropriate 
in different contexts. The following definition is standard for the controlla- 
bility of the composite system S + A. 

Definition 3.1 The system 5 + A is called completely controllable (CC) if, 
for any special unitary evolution P £ SU{nsnA)u there exists a time T > 
and a set of control functions, ui{t),. . . ,Um{t), defined in [0, T], such that 
the solution of ([7]) satisfies X{T) = P. 

In dealing with indirect controllability of S, it is important to take into 
account two processes: (i) the preparation of pA] (ii) the dynamical control 
of the full system S + A. The process of preparation of pA is typically an 
irreversible process completely separated from the unitary evolution of the 



* We assume linearity in the control as in dipole interaction, but this assumption is 
not crucial, and the treatment that will follow can be adapted to more general situations 
with minor modifications. 

^We shall use the notation u(n) (su(n)) to denote the Lie algebra of n x n skew- 
Hermitian matrices {n x n skew-Hermitian matrices with zero trace) while U{n) {SU{n)) 
denotes the corresponding Lie group of n x n unitary matrices (unitary matrices with 
determinant equal to 1). 



total system. This preparation is performed before dynamically controlling 
the system. Therefore, in the definition of indirect controllability we need 
to include the information on the set of possible initial states for the system 
A. This set is denoted by A^ in the following definitions. 

Definition 3.2 A state p'g is said to be reachable from ps given Aa, if there 
exists PA G Aa, and a set of controls ui(t), . . . , Um{t), defined in [0, T], such 
that the solution X = X{T) of ([7]) satisfies 

p's = Tta [x{T)ps paX\T)) . (11) 

If TZ is the set of unitary evolutions which can be obtained for the total 
system, the set of reachable states from ps given Aa is 

Asips) ■■= {Xps paX^pa G A^,X g n}. (12) 

In the previous definition, As{ps) depends on A^, and the choice of pA 
in Aa is part of our control strategy. Special cases are when Aa consists of a 
single state, that is the accessor admits only a specific initial configuration, 
or when it consists of states that have special features. The case of pure 
states is particularly important, both from an application viewpoint, and 
because of some results we shall prove in this paper (cf. proposition 15. 2p . 
Nonetheless, we will not a priori put any restriction on A^^. 

Definition 3.3 The system S is called fully indirectly controllable (FIC) 
given Aa if, for every initial state ps, the reachable set As{ps) is the set of all 
density matrices of 5. The system 5 is called unitary indirectly controllable 
(UIC) given Aa if, for any given initial state ps, the reachable set As{ps) is 
the set of all density matrices unitarily equivalent to ps- 

By definition, FIC implies UIC given the same set Aa- Moreover, if FIC 
(or UIC) holds given A^, it also holds given any A'^ such that A a ^ A^. 
It is a known result in quantum information theory that, if the system is 
CC and ua ^ n%, the system S is FIC given the set A^ consisting of one 
pure state for Ai^i We shall prove in Theorem O that if ns = 2, ua = 2 is 
sufficient to have this property. 

If UIC or FIC hold independently of the set Aa, we shall refer to these 
properties as strong. For instance, strong UIC means that, independently 
of the initial state pA, for all Xs G SU{ns) we are able to find a control 



"This is equivalent to Aa being the set of all pure states since under the CC condition 
we can transfer the state of A between two arbitrary pure states. So this statement 
is equivalent to pure FIC. It is proved that an arbitrary completely positive map can 
be performed on S (dilation theorem for completely positive maps [20]) and, according 
to [26], every possible state transformation can be realized through a completely positive 
map. 



strategy such that the solution of ^ is Xsps^s- Complete controllability 
of the total system implies the strong UIC notion since every transformation 
of the form Xs (8) 1 is, in particular, a unitary transformation, and gives the 
desired state transfer ps — )• Xsps^s, independently of the initial state pA- 
For this reason, most of the investigations on indirect controllability have 
focused on giving conditions for complete controllability (see, e.g., |15]). 
When UIC or FIC holds for the set A^ of pure states for A we shall refer to 
these properties as pure. Clearly strong implies pure for both properties. 



Remark 3.4 We remark that all the notions of indirect controllability which 
have been introduced, and the criteria which will be presented in the follow- 
ing, are invariant with respect to unitary changes of coordinates which affect 
the system S and A separately (local transformations) . They are not invari- 
ant under general unitary changes of coordinates that affect both systems 
S and A. 

According to classical results of geometric control theory |16] applied to 
quantum mechanics [8], [18], [21] . the set of all available unitary transfor- 
mations on the total system 5 + ^4 is dense in the connected Lie group e^ 
corresponding to the Lie algebra C generated by the set 

T := {iHs + iHi + iHA, il (g) Li, . . . , il (^ L^}, (13) 

where Hs^ Hj and Ha were defined in Q - ([6]), and it is equal to e^ if 
e^ is compact HI Therefore the Lie algebra C is of crucial importance in 
determining the structure of the reachable set for the total system and the 
system S. In particular the set As{ps) defined in ([12]) can be written as 

Mps) = {p'slpA G Aa, X G e^}. (14) 

The Lie algebra C is referred to as the dynamical Lie algebra associated 
with the system. We shall also consider the control Lie algebra, B which is 
generated by {il Li, . . . , il (5D L^}; it turns out that C is generated by B 
and iHs + iHj + iHA- 

For a subspace V of u{nsnA), let Tr^(V) denote the image of V under 
Tr^ in u{ns)- Consider an initial state for the total system, ps pA, and 
the subspace of u(nsnA), 



V := ad^ ispan{ips pa}) , (15) 



k=0 



^We shall not insist in the following on this distinction and in fact all the Lie groups 
we shall encounter are compact. 



where a,dj/P = V for any space V and, recursively, adj/P := [£, ad^ V]u 

Theorem 1 Consider the nontrivial case ps ^ — Ing- Assume that for all 
X G SU{ns) there exists C/ S e^ such that 

TvAiUps ® paU^) = XpsXl (16) 

ThenTvAiV) =u{ns). 

Proof. Since V is invariant under £, it is also invariant under the action of 
e^ by conjugation, i.e., Q E V =^ UQU^ G V for any [/ G e^. Now, from 
(fT6|) we have that, for every X G SU{ns), XipsX'^ G Tr^(V). This implies 
that 

SY>an{XipsX^\X G 5t/(n5)} C TrA(V). (17) 

Since ps / ;^lns, span{XipsXt|X G 5C/(ns')} = u{ns), it follows u(n5) = 
TV^(V). D 

Theorem [1] can be applied without calculating V when we recognize that 
ips®PA belongs to a subspace V invariant under commutation with £, which 
might in general include properly V. If Tr/i(V) 7^ u{ns), we can exclude the 
UIC property. One case is when ips ® pA ^ i^ ® span{zl„g ® In^} ^^'^ 
Tr^(£ ® spanjil^c; ® In^}) / ^{i^s)- An application of this will be used in 
Theorem m Another possibile case is when V = £-*" ® i'^ns ® Ia^a' which 
is also invariant under C. We have found an example showing that the 
converse of Theorem [1] is not true in general and it is an interesting open 
problem to understand how the condition of the theorem can be modified 
to have a necessary and sufficient condition. 

4 THE DYNAMICAL LIE ALGEBRA FOR THE 
INDIRECT CONTROL OF TWO QUBITS 

In this and the following section, we explore the case in which target S 
and accessor A are both two- level systems. We start with a study of the 
dynamical Lie algebra in all cases. We shall treat separately the two cases 
in which the control Lie algebra B is the full 1 ® su(2), and in which it 
is one-dimensional. These are the only possible cases since B is always a 
Lie subalgebra of (a Lie algebra isomorphic to) su(2) and there are no two- 
dimensional subalgebras in 5u(2). We shall also assume throughout that the 
interaction Hamiltonian Hi is non-zero since the case of trivial interaction is 
clearly non controllable and consists of the two systems evolving separately. 



*The space V can be calculated recursively from a basis of £, and because of u{nsnA) 
being finite dimensional, there exists a k such that V = ©t-^g ^^c (span{ipg ® Pa})- The 
main (in fact, defining) property of V is that it is the smallest subspace of uinsriA) which 
is invariant under £ and contains ips ® pA- 



4.1 Algebra with Pauli matrices 

We recall for future reference the definition of the Pauli matrices 



1 oh''y'-=(-i ' ^^ ^= -1 



-.-(r n),-.:=( . ;),-.:=(; .]■ (18) 



The related matrices 

o'x-=-o-x, ay:=-ay, (Tz:=-(Jz, (19) 

span the standard two dimensional representation of su(2) and satisfy the 
commutation relational 

[0'x,0-y] = az, [o'y,Crz] = CTx, [(Tz^CTx] = CFy, (20) 

and the anti-commutation relations 

Wj,crk} = -^Sjkl, {cTj, 1} = 2(Tj j,k = x,y,z. (21) 

Anti-commutation relations are useful when using the formula 

[A^B,C(^D] = ^ [A, C]0{B,D} + ^{A, C} [B, D] (22) 

in calculations. 

4.2 Model for two interacting qubits 

We specialize to the two qubits case the general model (JU, (0), (0), and 
write 1 = I2. The Hamiltonian of the system S has the form Hs = -B (8) 1, 
with iB £ su(2), and we can make a change of coordinates (cf. Remark 13. 4p 
to diagonalize B. This does not modify the form of Hj nor the dimension of 
B. The Hamiltonian of the accessor system alone has the form Ha := 10C 
with iC G su(2). Therefore we have 

iHs ■■= ojscTz ® 1, (23) 

iHi := i(Ta ®crx + i(^b (K) cr^, + iac (8) cTz, (24) 

iHA:=il®C, iCGSu(2). (25) 

Here and in the following, cTa, with a ^ x,y,z indicates a matrix determined 
by a vector a := [axjUy, azY S R '^ and given by 

aa--= ax(yx + ayay + azCTz, (26) 



^The commutator [-B, C] is defined as [B,C] := BC — CB and the anti-commutator 
{B, C} is defined as {B, C} := BC + CB. 



that is a general matrix in su(2)|^ Define the matrix 

K:=lh], (27) 




and write it as 

K:={D F), (28) 

with D of dimension 3x2 and F of dimension 3 x iP^I 

In this notation the model is determined by the matrix K, the parameter 
Us, the skew-Hermitian matrix iC and the control Lie algebra B. We are 
interested in the dynamical Lie algebra generated by the set 

S := {usCTz (g) 1 + ida (g) (Ta; + iaf, (g) ay + iac (Jz + il (^ C , B}. (29) 

In the following subsection we characterize C in all cases when B is the 
full Lie algebra 1 (8)su(2). Then, we shall give a necessary and sufficient 
controllability condition for the case where B is 1-dimensional and cos = 0, 
that is, the system S does not have dynamics by itself. 

4.3 The dynamical Lie algebra in the case B = su(2) 

Theorem 2 Assume dimB = 3. Let K, D, F and los defined as above. 
Then we have the following cases for the dynamical Lie algebra C generated 
by S in (2§(l . If cos ¥" 0; ''^^ have: 

la) If D ^ and F ^ 0, the system is completely controllable, that is, 
/: = su(4). 

lb; IfD^OandF = 0, then 

C = spanjcTz (g) 1, 1 1^ 5u{2), i{ax,cry) (gsu(2)}, (30) 

which is W-dimensional. 

Icj If D = and F y^O, then 

C = span{io-2 (g) su(2), ct^ (g 1, 1 (g su(2)}, (31) 

which is 7 -dimensional. 



^"in this correspondence between vectors in R "^ and elements in su(2), we have a^ -f-^ i, 
ay <r^ j, az ■'r^ k. This correspondence between su(2) and R "^ is, in fact, a Lie algebra 
isomorphism, where in R ^ the Lie bracket is replaced by the cross product. There is a 
correspondence also between inner products in su(2) and R ^ and orthogonal vectors in 
R '' correspond to orthogonal matrices in su(2). 

^^Notice that since we have assumed nontrivial interaction at least one between F and 
D must be nonzero. 
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If '^5 = 0, we have: 

2a) //rank(iC) = 1, there exists a £ su(2) such that 

£ := span{i(T(g)Su(2),l(g)Su(2)}. (32) 

2bJ If rank(i^) = 2, C is conjugate (and therefore isomorphic) to the Lie 
algebra C of point lb) above. 

2cj If rank{K) = 3, the system S + A is completely controllable. 

The condition rank(i(') = 3, or equivalently det{K) ^ 0, is the complete 
controhabihty sufficient condition of [15] for both the cases us = and 
cos ¥" 0- It is clear that there are cases included in la) where this condition 
is not verified and the system is nevertheless completely controllable. 

It is interesting to analyze the structure of the Lie algebra £ in the 
cases which are not completely controllable, in particular, in cases lb) and 
2b). This Lie algebra is isomorphic (and therefore conjugate, see Lemma 4.2 
in [3]) to sp(2) as it can be shown for example by using the test of Theorem 
8 in [3]. Therefore the total system 5 + A is pure state controllable, which 
means that every transfer is possible between pure states in S A. This 
implies that, if ps (d PA is a pure state, and the final desired state for ps 
is also a pure state, this transfer will be possible in the indirect control 
scheme. We shall see in Proposition 15.21 that a stronger property actually 
holds in this case, that is. Pure UIC. We shall consider this case also in the 
Example 15. H in a setting of interest for applications. We now give the proof 
of Theorem [2j 

Proof. We notice that since ;B = 1 (g) su(2) we can cancel 1 (gi C in ([25]) with 
an element of i3. We also notice that, in any case, the dynamical Lie algebra 
C is the same as the Lie algebra C generated by the set 

S' := {u>s(Tz'S)l,io-a'S>5u{2), 

tab (S> su(2), i<Tc (g) su(2), 1 (g) su(2)}. (33) 

The fact that £ C £' is obvious since the generators of £ can be obtained as 
linear combinations of the generators of £'. To show the converse inclusion 
we need to show that every generator of £' listed in (p3]) is in £. To do 
this take the Lie bracket of uso'z <8) 1 + iaa <^ cTx + ioh ® Oy + iuc ® az £ C 
with \®(Tz. This gives —iaa ®cry + ia^ ® Ox- Taking the Lie bracket of this 
with \® Gx gives ioa ® Oz and with Lie brackets with 1 ® ax,y and linear 
combinations, we obtain all elements of the form iaa ® <7, with arbitrary 
a E su(2). Analogously we obtain elements in iai, ® su(2) and iuc su(2). 
We therefore analyze the Lie algebra generated by the set in ()33p in all cases. 
Assume first ujs ^ 0. Assume both D and F are different from zero. 
Since D ^ 0, there exists one among aa, CTb and Uc with a nontrivial compo- 
nent perpendicular to Uz- Assume, without loss of generality, that it is Ua- 
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We can perform a rotation (i.e., a unitary similarity transformation) on the 
first qubit so that this component is proportional to (jy without modifying 
Gz, that is, we can assume without loss of generality that a a = otUy + /3a z, 
with a 7^ which we can therefore assume equal to 1. Taking the Lie bracket 
of iaa fX" o" for any a £ su(2) with a^ (E> 1, we obtain therefore iax cr. Taking 
the Lie bracket with a^ (8) 1 again, we obtain iay (8) cr. This means that the 
span of iax <8> a and iay (8> (T is in C Now since F / there is one of the aa^^c 
which has nonzero component along az (recall that our change of coordi- 
nates did not modify the components along cr^ since it amounts to a rotation 
in the x — y plane). Therefore we have every element of the form iai (gi a2, 
(Ti G su(2) and a2 € su(2). The remaining terms of the standard basis of 
5u(4) can be obtained by Lie brackets of the type [iai az,ia2 az], for 
appropriate ai and a2 in su(2). In the case D ^ but F = 0, the argument 
is the same but it stops at the point where we have generated all elements 
in the Lie algebra described in ([301) . If -D = but F / 0, there is at least 
one element among aa, a^ and ac which has nonzero component along az 
and this is the only possible nonzero component. Therefore we have already 
all the elements listed in (|31|) which span a proper Lie subalgebra of 5u(4). 
Consider now the case where uog = 0, i.e., there is no natural dynamics 
of the target system S. If rank([L> F]) = 1 then there exists a a G su(2) 
such that aa,b,c = ka,b,c(^, for coefficients ka,b,c- Therefore we obtain the Lie 
algebra described in ([32]) . If the rank is 2 there are two among aa, a^, and ac 
which are linearly independent while the third one is a linear combination 
of these two and therefore the corresponding matrices M ® a, with M free, 
do not contribute to the Lie algebra. Assume, without loss of generality 
that the linearly independent matrices are a a and o"ft. By a rotation again, 
up to a nonzero unimportant proportionality factor, we can assume that 
Ca = crx and (Tfe = cjy + aax- Therefore i(o"^, ay) ® su(2) is included in the 
Lie algebra C. By taking the Lie bracket \iax (8> az,iay ® az] one obtains the 
remaining basis element cr^ (8 1 in ([30p . In the case 2c) (rank 3) one adds 
to the previously generated elements another set of elements iac ® a with 
arbitrary a G su(2) and notice that ac must have, in the chosen coordinates, 
non zero component along az- Therefore one obtains as above the remaining 
elements in the standard basis of 5u(4). Alternatively, and more quickly, one 
can use the fact pointed out after the statement of the theorem, that the Lie 
algebra ([30]) is conjugate to 5p(2) and every such Lie algebra is known to be 
a maximal subalgebra [14J in su(4), that is, the addition of every nonzero 
element outside the Lie algebra but inside 5u(4), causes the generation of 
the whole Lie algebra su(4). D 
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4.4 Complete controllability when B is l-dimensional and 

ios = 

In the case where the control Lie algebra B is one-dimensional, that is, there 
is only one independent control, the analysis becomes more complicated. 
Among the other issues, one has to consider that the independent dynamics 
of the accessor system (i.e., the Hamiltonian Ha in dH)) cannot be canceled in 
general by the action of the control and play a significant role in determining 
the dynamics of the total system. Nevertheless, we can use the result of 
Theorem [2] jointly with the fact that, with the same parameters for the 
system, the dynamical Lie algebra in the case of B l-dimensional is always 
a subalgebra (not necessarily proper) of the corresponding dynamical Lie 
algebra for the case of B 3-dimensional. We restrict ourselves to the case 
where ujs = 0- 

Theorem 3 Assume los = and dimB = 1. Then the system S + A is 
completely controllable if and only if the following two conditions are verified: 

CI) det{K) / 0; 

C2) The components ofTvAiiHfHc) andTvs{HA) perpendicular to Tys{Hc) 
are not both zero. 

Notice that condition CI is the same as the condition 2c in the dim;B = 3 
case. To that, we have to add the generically satisfied condition C2 to have 
a necessary and sufficient condition in the case of B 1— dimensional. 
Proof, h is useful to make a change of coordinates on the system A, so as 
to make B = span{l ® cTz}- In these coordinates, the part of iC in (|29p 
which is parallel to cTz can be neglected because it is already contained in B. 
Therefore, we can consider only the part perpendicular to cTz which can be 
taken equal to oJA^yy, for a certain real parameter uja- In these coordinates, 
we need to redefine a a, CTb ^-nd ac and therefore K. Moreover we can make a 
change of coordinates on the system S to make ab = (ic^y, aa = {aa^ +^o'y) 
and (Tc = xax + ycry + zaz, for real coefficients a, /?, 7 and x,y, z. 
Conditions CI and C2, in these coordinates, become 

cr) 

/a 7 0\ 
det /3 7^ 0. (34) 

\x y z ) 

C2') 

wi + j;2 + y2 _^ 0. (35) 

The condition (j34p is necessary because it is needed in the B 3-dimensional 
case. If condition ([35]) is not satisfied, the set S in ((29|l can be taken, in 
these coordinates, 

S := {iaax (8) <t^ + ijCTy (^^ a^ + i(iay ® ay + iza^ cr^, 1 a^}, (36) 
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and both elements in this set are in the Lie subalgebra 

icTy (g) ay, tax's) crx,i(Jz'S^(Jz}- (37) 

We now prove the sufficiency of the conditions CI and C2 for complete 
controllability. We examine the Lie algebra generated by the set 5 which 
can be written in the given coordinates 

S := {iaa ax + iaj, (g) dy + iac (8) cr^ + uja'^ (^^ ay, 1(g) cr^}. (38) 



We notice that the Lie algebra generated by S in ()38p is the same as the Lie 
algebra generated by the set with three elements S' defined as 

S' := {iac (g o-z,i(^a 'g o^x + iab (8> cr^ + ujaI ^ay,l^ cr^}. (39) 

This can be easily seen because, if we take the Lie bracket of the first element 
of 5 in (j38p with the second element {l(giaz) two times, we obtain the second 
element in S' in (|39p . and by subtracting this from the first element of S, 
we obtain the first element in S' . We first characterize the Lie algebra C" 
generated by the second and third element in S' in ()39p . i.e., the Lie algebra 
generated by 

S" := {iaa (^ (^x + ic^b (^ (Ty + ujaI ® ay,l® az}. (40) 

Then we examine what happens when we include the first element of S' and 
see that, under the given assumptions, we generate the whole Lie algebra 

5U(4). 

In the given coordinates, we have 

S' := {i{xax+yay + zaz)'gaz,l'gaz, (41) 

i{aax + jay) (g ax + ifiay ® ay + ujaI ® cry}, 

and S" is given by the last two terms listed in ()4ip . 
Define 

k:=a^ + Au}\, (42) 

which is always non- vanishing since a ^ from ()34p . 
Define the following matrices: 

r^ := i iay ^ax± -j=(aiax ay - ujaI ^cTx)) , (43) 

T^ ■.= --(l^az± ^(«^2 «> 1 - ^ujAicTy ^az)j , (44) 



-.lb ._ 



ay =f —i={aiax ax + ujaI ® (7y) ] . (45) 
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It is straightforward to verify that 

[^ X 1 ^ y 1 -'-2) [^ y 1 ^ z 1 ^ X 1 [^ z 1 ^ X 1 y ^ \^^) 

and 

[r-,r^] = r;, [r^-,r;]=r-, [r;,r-] = r-. (47) 

Therefore, the Lie algebras C'^ := span{r+, r+, F^} and C~ := spanjF^:, F^^, F^}, 
are both isomorphic to su(2) (cf. ([2n|) ). and one can verify that 

[£+,£-] = 0, (48) 

so that the sum C^ © C~ is, in fact, a direct sum of Lie algebras. Consider 
the generators of £". We have 

Li = l®a, = -(F++F-), (49) 

and 

L2 = iaax ® (Jx + ilCTy (^ cr^, + i/Say (gi <Ty + uaI <8) (Ty 
^ [7(r+ + r;) + Vk{T- - F+) + /3(F+ + f; 



2L 



(50) 



From these expressions it follows that C" is always a subalgebra of C~^ © £ . 
Using (06]), (07]), dMD, we calculate 

+ f7' + (/3 + ^)')r-l, (51) 



and this, along with the fact that /3 and k are both different from zero, 
and the expression of Li in (jl5|) . shows that both F+ and F~ belong to C". 
Using this fact, it is easy to verify, using (06]) and (07]), that there are only 
three possibilities: 

1. 7 = and /3 = \/k and C" is given by 

£" = span{F~,Fy ,F^,F+} (52) 

2. 7 = and /3 = —\/k and C" is given by 

£" = span{F+,F+,F+,F;} (53) 

3. At least one between 7 and \P\ — yk is different from zero and C" is 
equal to £+ © £~. 
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The task is to examine, under the given assumptions, what happens when 
we add to these Lie algebras the matrix ixax 0(7 z + iydy (gi o"^ + izaz ®(Jz '■= 
iGc (8) (Tz ■ Let us assume without loss of generality that ac is scaled so that 
{(Tc,o-c} = —^1 (cf. ([2T]) ). The case ([52|) is considered in the Appendix. 
The case (15311 is similar. 



In the case where C" = £+ © £ , we calculate, defining T := ixax (8> <Tz + 
iycTy (g) az + izaz ® cTz, 

[(r+ + r;),r]=|i®a,. (54) 

If y 7^ 0, the result follows as before by comparison with the case of B 3- 
dimensional. If y = at least one between x and ua is different from zero. 
We calculate 

-^ [8uA[rt,T]+axiTt-T-),T] = 



Vk 



X + ^^' ) 1 ^ ^y, (55) 



and the result follows as before by comparison with the case of ;S 3-dimensional. 

n 



5 INDIRECT CONTROLLABILITY FOR THE 
TWO QUBIT SYSTEM 

The following theorem says that strong unitary indirect controllability is 
equivalent to complete controllability for the system of two qubits. Nonethe- 
less, we shall see that weaker notions of indirect controllability are still pos- 
sible without complete controllability. 

Theorem 4 A system of two qubits with B = l®5u(2) is strong UIC if and 
only if it is CC. 

Proof. We have already noticed in Section [3] how complete controllabil- 
ity implies strong UIC. Assume now that complete controllability is not 
verified. The dynamical Lie algebra is (modulo similarity transformations 
acting separately on the two qubits) one of the Lie algebras ([30]) . ([3T]) or 
([32]) . Therefore, it is sufficient to consider these cases. In the case ([3T]) . 
consider for S an initial density matrix of the form ps = '^ + Kiaz- Every 
element in C has the form 

A = 1^ ai + aaz ^1 + iaz a2, (56) 
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with a a real number, and o"i and o"2 elements of 5u(2). Since e is a compact 
Lie group, every element X can be written as the exponential of a matrix 
A £ C Therefore we have 

Xps paX^ = e^ps PAe'"^ = Z] ~ ^^aPs ^ PA- (57) 

Now, by induction on j, given the form of ps as a diagonal matrix, it is easily 
seen that ad-^^ps0pA is always the sum of a finite number of tensor products 
where the first factor is a diagonal element. Taking the partial trace with 
respect to A we still obtain a diagonal matrix. Therefore, from diagonal 
density matrices we can never reach matrices which are not diagonal. l£| 

In the case of C as in (j32p . the proof is the same. In fact, the Lie algebra 
(|32p becomes a subalgebra of (|3ip after a change of coordinates on S which 
diagonalizes a. 

In the case of C as in ([SH]) . consider ps = ^(1 + kiaz) for some real k. 
With this choice, ips®PA belongs to £ = (span(il(g)l))©£ which is invariant 
under C. The result follows by applying Theorem[Tl as Tr^(£) ^ u{ns) (see 
the discussion after the theorem) P^l D 

The proof of the Theorem in the cases (|3ip . (|32p works to show that 
UIC given a set A^ is not verified, even even if we fix the set Ka a priori. 
It is interesting to investigate if this is the case also for pop . In the next 
proposition we will show that if A^ is chosen as the set of all pure states, 
then the UIC property given A^ is verified in the case ([30]) . We first present 
an example in order to physically motivate the study of the Lie algebra ()30p . 
This example also illustrates the use of Lie group decomposition techniques 
in a direct analysis of the reachable set (fH 



Example 5.1 Consider the system of two qubits S and A interacting via 
an Ising interaction and with full control on the system A. It is described 
by 

iHs := 0-^ O 1, (58) 

iHc := Ux{t)l®ax + Uy{t)l® Gy, (59) 

Hi := ay ay. (60) 



^^Notice that this argument also holds if we consider specific sets Aa for pA- Therefore 
UIC given Aa is never verified, no matter what Aa is. 

^^In the case (jSOp . £ is conjugated to the symplectic Lie algebra sp(2) and therefore 
any transfer is possible between pure states. If we want to find an example to show that 
strong UIC is not verified we have to use ps ® Pa non-pure. 
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An application of Theorem [2] (or a direct calculation) shows that the dy- 
namical Lie algebra £ is given by (|30pPi To parametrize the elements of 
the corresponding Lie group e^, we use the Cartan decomposition ^| of £ 



with 



C = }CeV, (61) 



/C := span{(T2 (g) 1,1 (g) (T^jicJa; (gxT^jicTj^ (g) (T^}, 
V := span{l (g) cr^:, 1 (g ay, ia^ (g o-x,iax <g ay, 

iay® ax,iay®ay}, (62) 

satisfying 

[/C,/C] C /C, [V,V] C /C, [P,/C] C V. (63) 

Every element X in e can be written as X = K1AK2 where Ki and K2 
are in the Lie group e'^ and Jl is in the Lie group c^^ where ^ is a maximal 
Abelian subalgebra of C included in V. Since K. is isomorphic to u(2), we 
can use an Euler decomposition for Ki and K2 and write 

J^ — gtlO-z®lgt2l®0-Zg«t30"a;(g)(T2gt40-2(gll /g^\ 

_^ _ gSlO-z®lg«S20-a;®(TzgS30-z(g)lgS4l®0-2 /gg\ 

for real coefficients fi, . . . , ^4, si, . . . , S4. For the Abelian subalgebra A, we 
can choose 

^ := span{l (g cJa;, icXa; (g da;}, (66) 

so that 

i _ gail®a"a:gia20"3;(gl(Ta; /gyN 

for real coefficients ai and 02- Notice that we have used the freedom in 
ordering the factors in a different way in ()64p and ()65p . so as to minimize 
the number of parameters which are relevant for the indirect controllability 
problem. The first factor in ([M|) can be pulled out of the partial trace (cf. 
(|14p and ([3])) while the second factor is a unitary local operation on A which 
does not affect the partial trace and can be neglected. The last two terms 
in ([65]) can be included in the initial state ps ® PA- In conclusion, we have 
to consider a transformation ps ^ PA ^ Y Ps ® PA^^ with the matrix Y of 
the form 

Y = g«i30"a;®crzg*40-z(X)lgail®a-^gia2CTj;(g)a-2, Sio-z(X>lg«S2CTa;(J3(Tz /goN 



^■^When analyzing a quantum mechanical system whose dynamical Lie algebra £ is not 
the full su(n) one approach is to first decompose it according to the Levi decomposition [T2] 
into simple and Abelian Lie subalgebras [11]. In our case £■ is isomorphic to sp(2) and 
therefore it is simple. Levi decomposition only gives a simple component and no Abelian 
subalgebras. 
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The reachable set starting from ps ® pA will be 

^s{ps) = {e*i'^-(TrAfl)e*i'^^|ti,S3,54e R}, (69) 

Vt := (Ye'-"'^pse-'-''''®e'^''^pAe-'''''Y^'' 



with Y varying as in (j68p . Each term of Y in (j68p can be expressed in terms 
of a basis of m(4). We have for example 

e**^''"®''" = cosail(g)l + isinaia-^(g)a2, ai := -^, (70) 

and the other expressions are collected in the Appendix, see ()80p . Combining 
all the expressions for the factors in ()68p . we obtain for Y a formula of the 
type 

Y ■.= Co®l + C^(dd^ + Cy®dy + C;,®az, (71) 

for matrices Cq, Cx, Cy, and C2, functions of 6 real arbitrary parameters 
ai,. . . ,aQ. We have calculated the expressions of these matrices and the 
results are reported in the Appendix. 

Formulas such as (|69p. ()7ip allow in principle to calculate an expression 
of the reachable set for every initial state but in practice are typically very 
complicated and are useful only with the help of numerical simulations. We 
have plotted the reachable set starting from a state 

PS® PA= ( 2 (■*■ + '^^^^ + '^•^^^) ) "^ ( 2 *'"'" "^ "■''^^-' ) '^'^^^ 

in several situations, see Fig. (TJHl This illustrates the variety of cases that 
can arise depending on the underlying algebraic conditions. 

Proposition 5.2 The system of two qubits with Lie algebra C in ()30p is 
pure UIC. 

Proof. Inspection of the basis of C in ()30p tells us that the Lie group e 
contains elements of the form 

Z(^l, 10P, e'''^^''^\ (73) 

where Z is an arbitrary element of the form e*^^" with a real, P is any 
arbitrary element in SU{2) and t an arbitrary real number. We start with 
an initial state for ps (X" PA with pA pure. By applying a transformation 
10 P, we can always assume that 

P.= (J l).= E, (74) 
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Case 1 




Figure 1: Simulation of the reachable set as described in the main text, with Sx = 0, 
Sz — 1/2 and a^ = 1. The simulation consists of 729 points randomly distributed; the big 
spot marks the initial state. The state of the accessor system is pure and it is possible 
to reach points that are outside of the sphere corresponding containing the initial state. 
This means that it is possible to make the state more pure with the scheme studied here. 



Consider now a general matrix p'g eg) Ei and calculate 



e- a'^^p^ea'^- )^Ei= Xp'gX^ ® Ei, 



(75) 



where X denotes a generic transformation of the form e*^^" for a real. Using 
rS]) with p'g = ZipsZl, we have 



Tps®EiT^ = (^Z2XZipsZ\x^Z^^®Eu 



(76) 



and taking the partial trace, we obtain Z2XZipsZ\X'^ Z2- The claim follows 
because, from Euler decomposition of SU{2), every element of SU{2) can 
be written in the form Z2XZ1. D 

We remark that, given the structure of the Lie algebra C, we could 
have chosen a specific (arbitrary) pure state, and the proof of the previous 
proposition would have gone through. 
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Case 2 




Figure 2: Simulation of the reachable set as described in the main text, with Sx = 0, 
Sz = 1/2 and a^ — 0. The simulation consists of 729 points randomly distributed; the 
big spot marks the initial state. In this case the initial state is in a space invariant under 
the dynamical Lie algebra (modulo the identity it is contained in £) and remains in that 
space. This is reflected in the picture where only a line parallel to the z axis is flUed. 



5.1 Full Indirect Controllability 

We conclude our investigation of the indirect controllability of a system 
of two qubits studying Full Indirect Controllability. In general, complete 
controllability does not imply strong FIC since the state ps® PA = 3(1® 1) 
is a fixed state for the total unitary dynamics. If we consider the set A^ to 
be the set of pure states for A, we have however the following result. 

Theorem 5 Assume the system S + A is CC. Then S is pure FIC. 



Proof. We write the initial state of S as ps := Z],=i 2 ^ib)OI) with {\j)} 
an orthonormal basis in S, Yli=i2''^j ~ ^' ^^^ ^^^ initial state of A as 
PA '■= |^)(V'I- Therefore, 



PS® PA = ^rj\'4jj){il)j\, 



(77) 



where \iIjj) = \i) ® |^), and {ipj\ipk) = ^jk- An arbitrary unitary transfor- 
mation of the composite system, denoted by [/, will act as 

Ups®paU^ = Yl ^Mj)i^j\^ (78) 

i=i,2 
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Case 3 




Figure 3: Simulation of the reachable set as described in the main text, with Sx ~ 1/2, 
Sz = and a^ = 1. The simulation consists of 729 points randomly distributed; the big 
spot marks the initial state. 



with {(j)j\(f>k) = Sjk, and \(j)j) arbitrary orthonormal vectors vectors in S<^A. 
Now, consider the foUowing particular choices of U: 

1. [/ is a unitary operator such that every vector \(j)j) is maximally entan- 
gled. This is always possible, because there are 4 maximally orthonormal 
entangled vectors m S f^ A (cf. [2D]), and we need only 2 of them. In this 
case, 

Y,r,TrA{\cpj){^j\) = ll, (79) 

i=i,2 



p's 



that is, the initial ps is mapped to the maximally mixed state; 

2. U is the SWAP operator, that is Ups <?> paU = pa® Ps- In this case, the 

initial ps is mapped to the pure state |'i/')('0|. 

Since the unitary group is path-connected, there is certainly a continuous 
path in the space of states of the system 5, connecting the maximally mixed 
state to \tp){tp\, and representing reachable states. By further using local 
unitary operations acting on S", every final state can be obtained. D 
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Case 4 




Figure 4: Simulation of the reachable set as described in the main text, with Sx = 1/2, 
Sz = and Uz = 0. The simulation consists of 729 points randomly distributed; the big 
spot marks the initial state. The situation is similar to the one in Figure [5] since now the 
initial state is, modulo the identity, in i£^ which is also invariant under £. This explain 
why all reachable points are in the x — y plane. 



6 CONCLUSIONS 



In this paper, we have formally introduced the problem of indirect control- 
lability for quantum systems and we have given tools and tests to describe 
the reachable set. In particular, we have introduced physically motivated 
notions and proved the general Theorem [1] to conclude lack of indirect con- 
trollability. This theorem is a test at the Lie algebra level and therefore 
its application implies only operations from linear algebra. The case where 
both the target system and the accessor system are two level systems is the 
simplest non trivial case which is also of great importance in applications. 
To treat this system we have made full use of techniques in the Lie algebraic 
approach to quantum control. Our first step has been to describe the possi- 
ble dynamical Lie algebras that can arise (Theorems [2] and [3]). These results 
are of interest on their own, and generalize previous results in the literature. 
With this information at hand, we have proven in Theorem U] that we need 
complete controllability of the whole system in order to have unitary indi- 
rect controllability given any state of the accessor system. However, if we 
are allowed to consider only pure states for the probe, we can have unitary 
indirect controllability without complete controllability (Proposition 
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When the state of the accessor system is pure and we have complete con- 
trollabihty, we can have full controllability on the target system, i.e., move 
the state between arbitrary density matrices (Theorem [5]). 

This paper is a first step in the study of indirect controllability and sev- 
eral questions are of interest for future research. Among these, the extension 
of our results to systems with general dimensions and an analysis of how in- 
direct controllability of the system S through A depends on the relative 
dimensions of S and A. In this work we have given tools but also indicated 
the difficulties in pursuing this study including the fact that much of the 
criteria and results have to be necessarily dependent on the coordinates. We 
believe that a closer connection with the theory of entanglement in quantum 
information would be fruitful in further developing this important topic. 



7 Expression of Co, C^, Cy and Cz in example 15.1 

Recall that ai is defined in ([70]) . We have 

g<x,®it4 = C2I 1 + is2a, 1, aa := |, (80) 

^i^.m.a2 ^ C4I ® 1 + is4ax ® a^, 04:= -J, 

g-7,®isi ^ C5I ® 1 + issa^ ® 1, "5:=^, 
g*a.®..s2 ^ cf,l®l + ts(,a^®d,, "6:= -J, 

having defined Cj := cos(aj), Sj := sin(aj) for j = 3,4. We have for the 
matrices Co, Cx, Cy and Cz in example 15.11 

Co = C3C4C2+5C1+6I - S3S4C2-5C1+65-X 

-S3S4S2-5Cl-6^y + iC3C4S2+5Cl-6^z, 
Cx = iS3C4C2+5Cl-6l + iC3S4C2-5Cl-6ax 

iC3S4S2-5Cl+6ay - 830482+501+6^^, 
Cy = ZC3S4C2-5S1-61 + iS3C4C2+5Sl-6^x 

-i83C482+5Sl+6ay + C3S4S2-5S1+65-2, 
Cz = -1838402-581+6^ + iC3C4C2+5Sl+6ax 

-iC3C482+5Sl-6^y - S38482-5Sl-6az, (81) 

where we have defined Ci±j := cos(aj it qj), and Si±j := sin(aj it aj). 
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8 Proof of controllability in the case (1521) 



Using (j42p . we define cos{9) := -% and sin(0) := -^A. For brevity we denote 
by c and s, cos{9) and sin(6') respectively. The goal is to show that the Lie 
algebra generated by the matrices {r~,^^^,r+} and the matrix P := iuc^Oz 
is the full Lie algebra su(4). This is the Lie algebra generated by the matrices 

r~ := iay ^ ax - ciax ®ay + -l®ax, (82) 

- ^ 

T-, := iay ® ay + ciax ®ax + -1® ay, (83) 

Z:=l®az, (84) 

A := caz (8> 1 — 2siay (S" az, (85) 

P := iac az, (86) 

with the first four matrices forming a subalgebra. Applying the standard 
algorithm to calculate the Lie algebra generated by a set of matrices (see, 
e.g., [TO]) we first calculate all the Lie brackets of depth 1. We however do 
not report Lie brackets which give zero and Lie brackets among the first 
terms of (|82|) -(|85|) since they do not give any new directions. Moreover we 
also skip the Lie bracket with T~ . Since T~ is the Lie bracket of a linear 
combination of A and Z, say L, with F", from the Jacobi identity, for any 
matrix M we have [M,F-] = -[L, [Fj,M]] - [F", [M,L]]. So it is enough 
to consider Lie brackets with Z, A and F~ only. We calculate 

Qi := [P, A] = ic{yax ® az- xay ®az) + - {zax (8) 1 - xaz ® 1), (87) 

Q2 := 4[P, F^] = -yl (g) a^ + cxl ay - 2isac ® ax- (88) 

At step 2 we take the Lie brackets obtained at step 1 (Qi and Q2) with the 
generating matrices ([55j) - ([5U|) . We scale and eliminate directions that are 
already achieved. From [(5i,P], since z 7^ from CI', we obtain 

Ri=cac®l. (89) 

From [Qi,y4], we obtain 

R2 := ic zaz ®az^is yay (g) az - — (yaz ® 1 + zay (g) 1). (90) 

From [Qi,Fj], we obtain 

2 
c V sc ex 

R3 ■= ^ 1 (X" cjj^ - i—yax 'gax + — l'gax + 



+ i-{zaz'i^ay + xax'i^ay). (91) 
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From [(52, -P] we obtain 

-R4 := W<yc ®(yy + icxGc (g) cr^ + - 1 (g) (Ty. (92) 

From [Q2 ,Z], we obtain 

i?5 := yl CTy + cxl (gi (Ta; + 2isac (g Cy. (93) 

From [Q2>r~], we obtain 

Rq ■= -ic^xax <g 0-2 - iyay (g CTz + — {yaz g) 1 - zay (g 1). (94) 



Now calculate 



Si := [i?4, Z] = iyac ® (J^ — icxac g) o"y + - 1 g) ctj,.. (95) 



We have 



2scxySi - 2sy'^R4 + (c^x^y + y^)ii5 = (96) 

y^(y^ + c^a;^ — s^)l g) ctj^ + cxy{c^x'^ + y^ + s^)l g) cr^;. Now, under the given 
assumption, the second term in the above expression is equal to zero if and 
only if cxy = 0. In fact the coefficient in parentheses (c^x^ + y^ + s^) is 
not zero because this would imply s = 0, x = (since c^ = 1) and y = 0, 
contradicting assumption C2'. Therefore formula ()96p gives the needed extra 
local term unless at least one between c, x and y is zero. We only need to 
consider these cases. 

We notice that c ^ 0, otherwise CI' is violated. Moreover, if x = y = 0, 
we have z = 1, and s 7^ to fulfill C2'; in this case it is sufficient to consider 
[P, (52] = ~§1 ^ <^j/) which is the extra desired local term. It remains to 
consider the cases when x = and y 7^ 0, or viceversa. If x = 0, (|96p 
gives the needed extra local term unless y = is. In this case, the local 
term 1 (g a^ is obtained from [Q2,R4] by discarding the component along 
P. Finally, when y = and x 7^ 0, it is always possible to obtain the local 
term 1 ig) cjy by linearly combining Q2 and i?4, since c, s and x are all real 
numbers. 
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